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ABSTRACT 

We study the role of submillimetre galaxies (S MGs) in the galaxy formation process in 
the A Cold Dark Matter cosmology. We use the Baugh et al. semi-analytical model, which 
matches the observed SMG number counts and redshift distribution by assuming a top-heavy 
initial mass function (IMF) in bursts triggered by galaxy mergers. We build galaxy merger 
trees and follow the evolution and properties of SMGs and their descendants. Our primary 
sample of model SMGs consists of galaxies which had 850/im fluxes brighter than 5mJy at 
some redshift z > 1. Our model predicts that the present-day descendants of such SMGs 
cover a wide range of stellar masses ~ 10^° — IO^^/i^^A/q, with a median ^ IO^^/i^^A/q, 
and that more than 70% of these descendants are bulge-dominated. More than 50% of present 
day galaxies with stellar masses larger than 7 x lO^^/i^^Af© are predicted to be descendants 
of such SMGs. We find that although SMGs make an important contribution to the total star 
formation rate at z ~ 2, the final stellar mass produced in the submillimetre phase contributes 
only 0.2% of the total present-day stellar mass, and 2% of the stellar mass of SMG descen- 
dants, in stark contrast to the popular picture in which the SMG phase marks the production 
of the bulk of the mass of present day massive ellipticals. 
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1 INTRODUCTION 

Our understanding of the star formation history of the Universe has 
changed dramatically as new populations of star-forming galax- 
ies have been discovered at high redshifts. A key advance was 
the discovery of the cosmic far-infrared background by COBE 
dPuget etal.l (l99^. with an energy density similar to that in the 
UV/optical background, which implies that the bulk of star for- 
mation over the history of the Universe has been obscured by 
dust, with most of the radiation from young stars being repro- 
cessed by dust grains to far-infrared and sub-mm wavelengths. 
This was followed by the discovery of the submillimetre galax- 
ies (SMGs) using the Submillimetre Common User Bolometer Ar- 
ray (SC UBA) instrument on the James Clerk M axwell Telescope 
(JCMT) jSmail et al.lll997l ; iHughes et alj|l99^ , which are gen- 
erally interpreted as dust-enshrouded galaxies undergoing a star- 
burst. Follow-up studies have concentrated on SMGs with 850/im 
fluxes brighter than about 5mJy. Spectroscopy has revealed that 
the bulk of such SMGs lie at redshifts z ~ 1 — 3, with a median 
around z ^ 2, when t he Universe had only 20% of its current age 
jChapman et alj[2005l) . Explaining the abundance and redshifts of 
SMGs has posed a challenge to theoretical models of galaxy for- 
mation ever since their discovery. 

In the standard picture, due to the high star formation rates 
inferred in these galaxies (of ~ 100 — 10 00A/a?/ r~^ for a stan- 
dard initial mass function, e.g. lHughes et al.iil998i ; iChaDman et al.l 
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the SCUBA phase is expected to play an important role in 
the build-up of the stars in massive galaxies. Central to this topic are 
the questions: do SMGs represent the assembly of massive present 
day ellipticals? What is the typical duration of the SCUBA phase? 
What are the typical stellar and host halo masses of their descen- 
dants? Do the descendants all have high masses, and are they pref- 
erentially spheroid-dominated? And finally, what is the overall con- 
tribution of the SCUBA phase to the stellar mass of present-day 
galaxies? 

In this paper we address these questions using the predictions 
of a semi-analytical model for galaxy formation and evolution. A 
powerful feature of semi-analytical models is the ability to con- 
nect h igh redshift galaxies to their present-day descendants dBaugbl 
l2006l) . Since it is believed that the dust plays a central role in pro- 
ducing the submillimetre radiation observed from SMGs, the model 
must take into account the energy reprocessed by dust when the 
spectral energy distribution (SED) is calculated. 

The analysis in this paper is based on the lBaugh etai] (l2005l) 
version of the GALFORM semi-analytical model. iBaugh et alj found 
it necessary to assume a top-heavy initial mass function (IMF) for 
starbursts in order to reproduce the observed number counts and 
redshift distribution of the submillimetre galaxies, while also be- 
ing consistent with other observational constraints, including the 
luminosity function of Lyman-break galaxies at 2; ~ 3, the present- 
day optical and near- and far-infrared luminosity functions, and 
the gas fract ions and metallicities of present-day galaxies. In the 
iBaugh et al.l model, the SMGs seen in current surveys are mostly 
starbursts triggered by major and minor mergers between gas-rich 
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disk galaxies. This model was subsequently shown to predict evo- 
lution of g alaxies in the IR in good agreement with Spitzer ob- 
servations jLacev et a l. 2 008), and has also been used to make 
predictions for Herschel iLacev et alj|2010d) . The model also re- 
produces the local size-luminosity relation for lat e-type galaxies, 
and t h e fractions of early an d late-type galaxies jGonzalez et al.l 
I2OO9I) . [ Gonzalez et al. I ( l2010h perform a similar analysis to the 
one presented in this paper, but for Lyman-break galaxies, rather 
than SMGs. Lyman-break galaxies are identified using photomet- 
ric selection to isolate the Lyman-break feature present in the rest- 
frame ultraviolet spectral energy distribution of star-forming galax- 
ies (and which is accentua ted due to ab sorption by gas along the 
line of sight to the galaxy). iLacev et al.1 ( |2010b) present a detailed 
study of the nature of Lyman-break galaxies in hierarchical galaxy 
formation models. 

The semi-analytical galaxy formation approach has been used 
by ot her authors to model submillimetre galaxies. iKaviani et al.l 
mm were able to reproduce the number counts of SCUBA 
sources, but only by treating the dust temperature as a free pa- 
rameter, for which they needed to choose quite a low value (20- 
25 K), in appare nt conflict with subsequent observations of sub- 
mm SED shapes dCoDpin et al.l2008l) . lFontanot et alj bOOTh repro- 
duced the number counts of galaxies in the band using a 
self-consistent calculation of the dust temperature and a standard 
Salpeter initial mass function for all types of star formation, but 
their model does not match the observ ed redshift distribu tion of 
SMGs. lSwinbank et al.l ( I2OO8I) tested the lSaugh et akl j2005h model 
by comparing in detail the SEDs and stellar, dynamical, gas and 
halo masses of submillimetre galaxi es against observa tional data, 
finding broad agreement. Recently Dave et al. ( 2010) have pro- 
posed an alternative model of SMGs based on hydrodynamical sim- 
ulations, in which SMGs are massive galaxies undergoing quies- 
cent star formation with a normal IMF. However, they do not cal- 
culate the sub-mm luminosities of their simulated galaxies, and so 
do not show that their model can reproduce the observed number 
counts and redshift distribution of SMGs. Prior to the present paper, 
no attempt has been made to relate the SMGs predicted in models 
with their present-day descendants. 

This paper is laid out as follows. In Section 2, we outline the 
galaxy formation model GALFORM used to predict the properties 
of SMGs. In Section 3, we present examples from the model of 
galaxy evolution histories which produce SMGs. In Section 4, we 
study properties predicted for SMGs at different redshifts, such as 
stellar and host halo masses and morphologies, and also the du- 
ration of the SMG phase. In Section 5, we show the model pre- 
dictions for the present-day descendants of SMGs, including their 
stellar and host halo masses and morphologies. In Section 6, we 
compute what fraction of present-day galaxies had SMG progeni- 
tors, and in Section 7 we calculate the contribution of SMGs to the 
present-day stellar mass density. Finally we present our conclusions 
in Section 8. 



2 GALAXY FORMATION MODEL 
2.1 Basic components 

We use the Durham semi-analytical gal axy format i on m odel, 
GALFORM, which is described in detail by ICole et al.l ( I2OO0I) and 
iBenson et sA\ j2003h . The model is set in the Cold Dark Matter 
(CDM) cosmology, with density parameter, f^o = 0.3, a cosmolog- 
ical constant, Q,\ = 0.7, Hubble constant -ffo = 70kms~^Mpc~^, 



baryon density, Q.b = 0.04 and a power spectrum normalization 
given by erg = 0.93. Galaxies in the model are assumed to form in- 
side dark matter haloes, and their subsequent evolution is controlled 
by the merging histories of the host dark matter haloes. The physi- 
cal processes modelled include: i) the hierarchical assembly of dark 
matter haloes; ii) the shock heating and virialization of gas inside 
the gravitational potential wells of dark matter haloes; iii) the radia- 
tive cooling and collapse of the gas to form a galactic disk; iv) star 
formation in the cool gas; v) the heating and expulsion of cold gas 
through feedback processes such as stellar winds and supernovae; 
vi) chemical evolution of the gas and stars; vii) mergers between 
galaxies within a common dark halo as a result of dynamical fric- 
tion; viii) the evolution of the stellar populations using population 
synthesis models; ix) the extinction and reprocessing of starlight by 
dust. 

The merger histories of dark matter haloes are calculated us- 
ing a Monte Carlo techniq ue, following the formali s m of extended 
Press & Schechter the ory dPress & Schechtei|[T974l ; lLacev & Cold 
ll993l ; IColeetalj2000l) . The formation and evolution of a represen- 
tative sample of dark matter haloes is modelled. 

We use the s ame galaxy formation parameters adopted by 
iBaugh et al.l ( |2005|), but we use a mod ified version of the Monte 
Carlo merger trees jParkinson et al.l2008il . which better reproduces 
the properties of dark matter halo merger tr ees extracted from the 
Millennium Simulation te pringel et alj2005 ). This is necessary be- 
cause we wish to connect high redshift galaxies with their present- 
day descendants, and so we need a prescription for building merger 
histories of dar k matter haloes w hich is accurate over a long inter- 
val in time. The lParkinson et al] algorithm for bu ilding halo merger 
trees is a modified version of that introduced bv lCole et alj 1 I2OOOI) 
and used in Baugh et al. (2005), and has been tuned to match the 
merger histories extracted from the Millennium Simulation. This 
modificati on does n ot significantly affect the published predictions 
of the Bau gh et"al] model, since those did not de pend on the a ccu- 
racy of the trees over long time intervals. Instead, feaugh et al] laid 
down grids of halos at a range of redshifts to compute the number 
counts of SMGs and the cosmic star formation history. We have 
checked that using the lParki nson et al.' tree algorithm does not alter 
the predi ctions presented in Baugh et al. 

The iBaugh et alj ^2005 ) model uses feedback from a super- 
wind to suppress the formation of massive galaxies and uses a 
top-heavy stellar IMF in starbursts triggered by galaxy mergers 
to reproduce the observed numbers of submillimetre and Lyman- 
break galaxies. An alternative model was presented by Bower et al] 
(2006), which adopts a standard IMF throughout, but uses AGN 
feedback to reduce the number of massive galaxies. However, the 
iBower et al.l model underpredicts the number counts of submillime- 
tre by more than an order of magnitude, so we do not consider it 
further here. Fo r a complete list of th e differences between these 
two models, see iGonzalez et al] ( l2009h . 

2.2 Mergers and star formation bursts 

When dark matter haloes merge, the galaxy in the most massive 
halo is assumed to become the central galaxy in the new halo while 
the other galaxies become satellites (note that the central galaxy 
identified in this way is not necessarily the most massive galaxy 
in the new halo). The model assumes that both central and satel- 
lite galaxies can form stars from their cold gas reservoirs, but only 
central galaxies can accrete more cold gas by cooling from the sur- 
rounding halo. dFont et all l20og) considered an alternative model 
in which the hot gas halos of satellite galaxies are only partly 
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stripped, but this has little effect on SMGs in our model, which 
are typically central galaxies.) 

Mergers in the model are assumed to happen only between 
satellite and central galaxies within the same halo. Bursts of star 
formation in the model are triggered by galaxy mergers. Two types 
of mergers are defined, major mergers and minor mergers, accord- 
ing to whether or not the ratio of the mass of the smaller galaxy 
to the larger galaxy M2/M1 exceeds the threshold /dup. Bursts 
are assumed to be triggered in all major mergers, but also in mi- 
nor mergers that satisfy M2/M1 > /burst and where the ratio of 
gas to stellar plus gas mass of the larger galaxy exceeds /gas.crit. 
The values adopted in the iBaugh et alj model are /diip = 0.3, 
/burst = 0.05 and /gas.crit = 0.75. Mergers can change the mor- 
phologies of galaxies. In a major merger, the stellar disks of the 
merging galaxies are assumed to be transformed into a new stel- 
lar spheroid (resulting in a pure bulge galaxy) and any gas present 
in the disks is converted into stars in the burst, or ejected by su- 
pernova feedback associated with the burst. In minor mergers, the 
stellar disk of the larger galaxy is preserved and the stars from the 
smaller galaxy are added to its bulge, so the merger remnant is a 
disk-l-bulge galaxy. If the condition for a burst in this minor merger 
is satisfied, then the gas of the smaller galaxy is converted into stars 
and added to the bulge, while if the condition is not satisfied, any 
gas is added to the main gas disk . 

As shown in 'Baugh et al] j2005h , our model predicts that 
most SMGs seen in current sub-mm surveys are starbursts trig- 
ger ed by major and minor mergers involving gas-rich disk galax- 
ies. iDekeletalJjiooa , and references therein) have proposed a pic- 
ture in which the most actively star forming galaxies at high red- 
shifts are fed by cold streams of gas falling into galaxy halos. Al- 
though this picture appears very different from our model, the dif- 
ferences are actually less dramatic than they appear at first sight. 
Firstly, the GALFORM semi-analytical model has always included 
a "cold" mode of accretion, in which gas falls to the halo centre at 
the free-fall rate, for halos in which the cooling time is less than the 
free-fall time for gas shocked at the virial radius cCole et a l. 2000). 
This " cold" mode dominates a t lower halo masses and higher red- 
shifts. Benson & Bower ('2010') have recently shown that using the 
iBim boim & Dekel (2003) criterion to distinguish "cold" and "hot" 
accretion modes instead of the standard GALFORM criterion slightly 
shifts the halo mass at which the transition from "cold" to "hot" ac- 
cretion occurs, but leaves gas accretion rates and SFRs in the model 
almost unchanged. They also find that the importance of the "cold" 
mode is gre atly reduced when realistic feedback is included. Sec- 
ondly, in the lDekel et al J picture, the bright SMGs are actually trig- 
gered by the infall of large clumps of cold gas rather than smooth 
streams - effectively gas-rich galaxy mergers. 

2.3 Top-heavy IMF in bursts 

As m entioned in the introduction, an important feature of the 
iBau gh et al. model is that stars are assumed to form with differ- 
ent initial mass functions (IMFs) in different environments. This 
differentiation depends on whether they form quiescently in disks 
or in bursts following a galaxy merger. In the case of disks, the 
model adopts a standard solar neighbourhood IMF proposed by 
[Kennicutt ( 1983), for which the number of stars produced varies 
with mass as dA''/dlnm oc m~^, with x = 0.4 for m < IMq 
and X = 1.5 for m > IMq. In the case of bursts, a top-heavy 
IMF is used, where s = 0. In either case, the IMF covers a mass 
range 0.15 < m < 12OAf0. This top-heavy IMF in bursts has two 
effects in the model: i) it increases the total UV radiation from the 



young stars formed and ii) it increases the metal production and 
hence also the amount of dust. 

When stars die, they return gas and metals to the interstellar 
medium (ISM). We treat this process using the instantaneous recy- 
cling approximation, in which the metal ejection is quantified by 
the yield, p, (which is the mass of metals returned to the ISM per 
unit mass of stars formed) and the gas return is specified by the re- 
cycled fraction, R, (defined as the mass of gas returned to the ISM 
per unit mass of stars formed). We have calculated the values of p 
and R for o ur two IM Fs based on the results of stellar evolution 
calculations jCole et a l. 2000; B augh et aklboOSh , and find values 
Pquicsccnt = 0.023, -Rquicsccnt = 0.41 for quiesccnt star formation 
andpburst ~ 0.15, -Rburst ~ 0.91 for bursts. 

There is independent support for the idea of vari ations in the 
IMF from a variety of observations, as discussed in 'Lac ev et"ail 
12008). We mention a few of these here. Fardal et al. (2007!) found 
a discrepancy between the observed ratio of the total extragalac- 
tic background radiation to the observed stellar density today and 
the value predicted using the observationally inferred star forma- 
tion history, if a standard IMF is assumed. They find that a more 
top-heavy IMF is needed (increasing the proportion of stars formed 
with 1.5 A/0 < m < 4M0) in order to reconcile these different 
observations. They argue that alternative solutions in which non- 
stellar sources make large contributions to the background light 
appear unlikely to resolve the discrepancy. Another s uggestion of 
variations in the IMF comes from the work of lvan Dokk um (200^, 
who compares the luminosity evolution of massive cluster galaxies 
(at 0.02 ^ 2 ^ 0.83) with their colour evolution. He finds that 
the evolution of the rest-frame U — V colour is not consistent with 
the previously determined evolution of the rest-frame M/ Lb ra- 
tio for a standard IMF, but shows that using a different IMF with 
a flatter slope compared with a standard IMF can help to solve the 
discrepancy. 

2.4 Reprocessing of stellar radiation by dust 

The IR/sub-mm emission detected from SMGs is believed to origi- 
nate from the dust heated by radiation from young stars. Therefore 
a treatment of absorption and reemission o f radiat ion by dust is a 
key element of the model. In iBaugh et al.l ( l2005h . we calculated 
the reproc essing of stellar radiation by du st using the GRASIL 
code (Silv a et al. I ll998l : iGranato et al.ll200(]|) . which calculates the 
distribution of dust grain temperatures within each galaxy based 
on a radiative transfer calculation and a detailed grain model. The 
GRAS I L model has been shown to be able to very successfully ex- 
plain the SEDs of spiral, irregular, elliptical and starburst galaxies, 
including ULIR GS, over a huge range of wavelength, from the far- 
UV to the radio jSilva et alj[l993 : IVega et alj|2008l: ISchurer et al.1 
[2009.) . However, a drawback of the GRASIL code is that it takes 
several minutes of CPU to run a single galaxy, and so it is very time- 
consuming to run it on large samples of galaxies. For the present 
paper, we need to follow the evolution (including the sub-mm emis- 
sion) of large samples of galaxies over many timesteps, in order to 
be able select the small fraction of galaxies undergoing an SMG 
phase at any cosmic time, and it is not computationally feasible to 
do this using the GRASIL code directly. 

We have therefore developed a simplified model of dust emis- 
sion which approximately reproduces the results from the GRAS I L 
code at long wavelengths. At the sub-mm wavelengths of interest 
here, the SED of dust emission is insensitive to details of the dust 
model such as the dust grain size distribution and the effects of 
small grains and PAH molecules, allowing us to apply a simpler 
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Figure 1. Left panels: galaxy merger trees for different z = galaxies. The size of the symbol is proportional to the stellar mass as indicated by the key. The 
red circles indicate galaxies in a bright SMG phase (when the 850/xm, flux exceeds 5 mJy). Right panels: evolution with redshift of the galaxy on the main 
progenitor branch. Top subpanels: evolution of the bulge-to-total stellar mass ratio (B/T). The dashed vertical line indicates the redshift at which the central 
galaxy (black line) becomes a satellite galaxy (magenta line) after a halo merger Lower subpanels: the symbols show the evolution of the stellar mass (blue), 
cold gas mass (green), dark matter host halo mass (black) and the mass formed in ongoing bursts (red) in units shown on the left axis. The red hne shows the 
evolution of the SbO/im flux, in units shown on the right axis, and the dashed vertical black lines indicate the redshift at which the galaxy undergoes a merger 
with another galaxy (see left panel). From top to bottom, the stellar masses of the galaxies at z = are Af, = 2.7 X 10^'', 1.0 X 10^^ , 2.9 X 10^^ h^^ Mq 
respectively. In the case of the upper two galaxy trees, we show all of the progenitors with stellar mass in excess of 10^h~^ Af©; in the lower tree, we plot 
progenitors with stellar mass greater than 10^h~^ Mq. 
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approach. As in GRAS I L, we assume that the dust is in two phases, 
molecular clouds and a diffuse medium, with stars forming inside 
the molecular clouds and then escaping into the diffuse medium on 
a timescale of a few Myr. We use an approximate radiative transfer 
calculation (different from the one in GRAS I L itself) to compute 
how much stellar radiation is absorbed in each dust component. 
The most important difference with GRAS I L is that in our sim- 
plified model we assume the dust temperature within each of the 
dust phases is uniform (so there are only two dust temperatures in a 
galaxy). In contrast, in GRAS I L each size and composition of dust 
grain has its own temperature, which depends on position within 
the galaxy (which determines the intensity of the stellar radiation 
field which heats it, via the radiative transfer calculation). Further- 
more, GRAS I L includes the effect of fluctuating grain temperatures 
for small grains and PAH molecules. Thus in a full GRAS I L cal- 
culation, there is a whole distribution of dust temperatures within 
each galaxy. A final approximation we make in our simplified dust 
model is to represent the dust emissivity by a power law in wave- 
length. These differences between our two-temperature model and 
GRAS I L are crucial when calculating the dust emission at mid-IR 
wavelengths (which is dominated by emission from small grains 
and PAH molecules, and by warm dust in star-forming clouds), 
but have much less effect at the sub-mm wavelengths of interest in 
this paper. We give more details about our simplified dust model in 
iLacev et alj feOlOa'), where we also show some comparisons with 
GRAS I L. In summary this model works remarkably well at long 
wavelengths in the Rayleigh- Jeans tail of the modified black body 
spectrum of the dust emission. 



3 GALAXY MERGER TREES : EXAMPLES WITH SMG 
PROGENITORS 

To follow galaxy evolution in GALFORM we output the model 
predictions at different redshifts. We label a galaxy as a bright 
SMG when its flux exceeds 5 mjy. As examples of dif- 

ferent merger histories which can produce bright SMGs, we plot 
in the left panels of Fig. [T] galaxy merger trees for three different 
present-day galaxies, covering a range in present stellar mass from 
~ 3 X lO"^/i"^A/0 to ~ 3 X lO"/!."^A/0. In these examples we 
require that a SMG phase appears in the main progenitor branch 
(leftmost branch), which we define as follows: starting at the final 
time, we step back in time (up the tree), following the branch of 
the progenitor with larger stellar mass each time there is a merger. 
(Note that the main progenitor need not be the most massive pro- 
genitor at all times.) We indicate with a red circle when a galaxy at 
a given redshift is a bright SMG and with a blue circle otherwise. 
In these selected examples, the more massive present-day galaxies 
can show more than one bright SMG episode. 

In the right panels of Fig. [T] we show the evolution of the 
850/im flux for the main progenitor galaxy (leftmost branch) in the 
trees. Following each merger that triggers a burst, there is a rapid 
increase in the 850/im flux, followed by a slow decline. The time 
interval when the flux is above 5 mJy is different in the different ex- 
amples. In the same panels, we also plot the evolution of the stellar, 
gas and dark matter halo masses for the main progenitor, together 
with the mass formed in ongoing bursts. In the upper sub-panels, 
we plot the evolution of the bulge-to-total stellar mass ratio B /T as, 
an indication of the morphology of the galaxy, and we also indicate 
when the central galaxy became a satellite galaxy. 

In the first example shown in Fig. [T] the central galaxy be- 
comes a satellite galaxy ai z = 2.3 (note the large increase in the 



host halo mass due to the merger with a more massive halo). In the 
second example, the galaxy becomes a satellite at z = 0.3 and in 
the third example, the galaxy becomes a satellite at z — 0.1. As 
mentioned previously, in the standard GALFORM scheme, satellite 
galaxies cannot accrete any more gas by cooling, so star formation 
can only continue until the pre-existing cold gas reservoir is used 
up, which limits the increase in stellar mass once a galaxy becomes 
a satellite (green line). 

We also indicate in the right-hand panels when the galaxy un- 
dergoes a merger with another galaxy. If there is a burst (see Sec- 
tion 2.3) the gas from the merging pair is used up to produce a 
burst of star formation (which raises the 850/im flux). Note that 
the galaxy mergers at 2: = 2.7 and z = 2.5 in the first and second 
examples, respectively, are minor mergers since the morphology of 
the resulting galaxy is a mix of bulge-l-disk (B /T ~ 0.5 in both 
cases). The galaxy mergers at z = 5.9 and z = 2.1 in the first and 
second examples, respectively, are major mergers since the result- 
ing galaxy is a pure bulge galaxy with B/T — 1. Since more hot 
gas can be accreted by cooling onto central galaxies, more star for- 
mation takes place, changing again the morphology of the galaxy. 
In the third example, the main progenitor becomes bulge dominated 
following a major merger at z ~ 1.5 and remains bulge dominated, 
even though it is still the central galaxy in its halo until 2 = 0.1. 
In this example, there are a number of mergers with other galaxies 
which add mass to the spheroid (indicated by the vertical dashed 
lines in the main panel on the bottom right). 



4 PROPERTIES OF SMGS AT DIFFERENT REDSHIFTS 

Having identified all the SMGs in the model, we can study the 
evolution of their properties. In the model , the redshift distribu- 
tion of SMGs shows a peak at z ~ 2 iBaugh et all see also 
Swin bank et al.l) similar to that observed by Chapman et alj j2003l . 
[2005). For convenience in what follows, we define two types of 
submillimetre galaxies according to a threshold in the 850/im flux. 
We define a galaxy as a bright SMG when its flux exceeds 5mJy 
and as a faint SMG when its flux exceeds ImJy. The faint sam- 
ple therefore includes the bright sample, but these bright galaxies 
only account for a small fraction of the galaxies in the faint sam- 
ple. Most currently observed SMGs have fluxes around 5mJy or 
brighter (and so are bright SMGs by the definition used here). We 
present the faint SMGs as a comparison sample. 

Our samples of SMGs are effectively volume limited. We con- 
sider dark matter halo merger histories for a grid of halo masses, 
generating many different examples for each mass on the grid. 
SMGs are identified according to their flux as outlined above, given 
the redshift of a particular branch in the merger history. The distri- 
butions of properties from different mass trees are combined by 
assigning a weight to each tree based on the abundance of haloes 
of that mass and the number of example histories produced. 

4.1 Duration of the SMG phase 

We define the duration of the bright SMG phase to be the time 
for which the 850/im flux of a galaxy is above 5 mJy. This time 
is typically smaller than the duration of the star formation burst 
which triggered the flux increase (see the examples in Fig.[T](. In 
Fig. |2] we plot the distribution of SMG phase durations all of the 
bright SMG phases identified in the redshift range 2 < z < 3. 
The distribution is that for a volume-limited sample. The durations 
range from < 0.04 Gyr to roughly 0.5 Gyr. The median of the 
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broad-band fluxes. However, this comparison is complicated by the 
fact that the model includes a top-heavy IMF in starbursts, while 
the observational estimates assume a solar neighbourhoo d IMF for 



all stars (see also the discussion inlLacev et al. 



recent observational analysis bv iHainline et al 



0.3 0.4 0.5 0.6 0.7 
duration [Gyr] 

Figure 2. Duration of the bright SMG phase, i.e. the time for which the 
flux is Sv{9,^0^lm) > 5m Jy, for SMGs at 2 < z < 3. The distribution is 
that for a volume-limited sample. The median of the distribution is 0. 11 Gyr 
(indicated by the arrow). The distribution is normaUzed to unit area. 



distribution is 0.11 Gyr. For comparison, the median duration of 
the star formation bursts that triggered the bright SMG phase is 
0.66 Gyr. These predicted durations are similar to obs ervational 
estimates base d on UV spectral fittin g dSmail et al.ll2003h . duty cy- 
cle arguments I Chapman et alj200^ and gas depeletion timescales 
jTacconi et alj200a) . 



4.2 Stellar and Halo Masses 

The left-hand panels of Fig. |3] show the median stellar and halo 
masses of SMGs identified at different redshifts (the right-hand 
panels of this plot show the properties of the descendant galaxies 
and are discussed in the next section). The medians are calculated 
for volume-limited samples of SMGs at each redshift. The top left 
panel of Fig. |3] shows the median stellar mass and host halo mass 
for the bright SMGs as a function of redshift. For the bright SMG 
sample at redshift z=1.5, the stellar and parent halo mass are respec- 
tively 4 and 2 times more massive than the bright SMGs identified 
at redshift z=6. For comparison, we show also the results for faint 
SMGs (i.e. those with 5*1^ > ImJy) in the bottom left panel. Faint 
SMGs identified at redshift z=6 have typically 1/4 of the stellar 
mass of the bright SMGs identified at the same redshift, and they 
live in haloes with 1/3 the mass. By redshift z=1.5, this difference 
in host masses of faint relative to bright SMGs has shrunk to around 
a factor 0.7 for both stellar and halo masses. 

At redshift z = 2, close to the median redshift of bright 
(Sv > 5mJy) SMGs (both in the model and in the observational 
data), the model predicts median stellar masses of 2.1 X 10^'' and 
1.7 X lO^°ft"^M0 for the bright and faint SMGs, and host halo 
masses of 2.2 x 10^^ and 1.5 x lO^^/i^^M© respectively. The 
predicted stellar and halo masses of S MGs in this model were com- 
pared with observational constraints in ISwinbank et ^ ( l2008h . The 
main observational constraint on halo masses comes from mea- 
surements of galaxy clustering. Swinbank et al. found good agree- 
ment between preliminary predictions of SMG clustering from the 
model and obs ervational data , but th is issue will be e xamined in 
more detail in lAlmeida et al.l ( l2010h . I Swinbank et alj also found 
that the predicted stellar masses of SMGs were lower than ob- 
servational estimates based on fitting stellar population models to 



dated stellar population models and correcting for AGN contami- 
nation of the rest-frame near-IR light, finds a median stellar mass 
of ~ 5 x lO^°/i"^M0, for a sample of SMGs with S„ ~ 5m Jy 
and 2 ~ 2, assuming a Kroupa (20 0l|) IM F. (Note that the typi- 
cal stellar masses found by Hainline et al.' are significantly lower 
than the values estimated by Michalowski et al. (2010) for a sim- 
ilar SMG sample, but not allowing for AGN contaminati on of the 
broad-band SED.) The median stellar mass found bv Hainline et al.l 
is not very different from the median value predicted by our model 
for bright SMGs sX z = 2, even before allowing for the difference 
in IMFs and star formation histories between the model and what 
is assumed in the observational analyses. We will investigate this 
issue in more detail in a future paper. 

4.3 Morphology 

In a similar way, we can study the morphologies of SMGs iden- 
tified at different redshifts. As mentioned above, a major merger 
produces a pure bulge galaxy, whereas a minor merger can result 
in a galaxy with a disk and a bulge (if the primary galaxy had a 
disk at the moment of merging). Note that in our model, the star- 
burst responsible for a galaxy appearing as an SMG starts after the 
two galaxies have merged. The left-hand panels of Fig.|4]show the 
model predictions for the B /T ratio of SMGs identified at different 
redshifts, and the right-hand panels, which we discuss in the next 
section, show B /T for SMG descendants. As in Fig. [3] the me- 
dians are calculated for volume-limited samples of SMGs at each 
redshift. In the top left panel of Fig. [4] we plot the median of the 
bulge-to-total stellar mass ratio B /T (where B/T=0 indicates a 
pure disk galaxy and B /T=\ indicates a pure bulge galaxy) for the 
bright SMGs. At high redshift, bright SMGs have intermediate val- 
ues of B/T {B/T ~ 0.4), but the typical B/T increases to ~ 0.7 
(indicating a bulge-dominated galaxy) as the redshift decreases to 
2 ~ 1, before dropping again to B/T ~ 0.4 at lower redshifts. At 
2; = 2, bright SMGs are predicted to have a median B/T = 0.6, 
indicating a midly bulge-dominated galaxy, but with quite a broad 
distribution, from 0.4 to 1 (10-90% range). 

For comparison, we show also the results for faint SMGs se- 
lected with > ImJy in the bottom left panel of Fig. |4] These 
faint SMGs typically have intermediate values of B/T at all red- 
shifts, with a very wide range (from 0.1 to 1), and the samples at 
higher redshifts have slightly higher values of B/T compared to 
those at lower redshift (showing an opposite trend as compared to 
brig ht SMGs). 

ISwinbank et alj JioTol) have rececently investigated the mor- 
phologies of bright SMGs at 2 ~ 2 using HST optical and near-IR 
imaging. They find that the light profiles are fit by Sersic indices 
n ~ 2 — 2.5, compared to the values n — 1 and n = 4 expected 
for exponential disks or r^'^^-law spheroids respectively. This re- 
sult indicates that these SMGs are mildy bulge-dominated, entirely 
consistent with the B/T values predicted by our model for SMGs 
at similar fluxes and redshifts. 



4.4 The SMG phase trigger: minor or major mergers? 

As we discussed before, most (> 99%) of bright SMGs in our 
model are starbursts triggered by galaxy mergers. Bursts happen 
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Figure 3. Left panels: The median stellar and host halo masses for SMGs selected at different redshifts. Right panels: median stellar' and host halo masses of 
the z=0 descendants of SMGs identified at different redshifts. Top panels: bright SMGs selected with Si, > 5mjy. Bottom panels: faint SMGs selected with 
Si, > Imjy. Filled squares show stellar masses and open squares show halo masses. EtTorbars indicate the 10% and 90% percentiles. The medians are all 
calculated for volume-limited samples of SMGs at each redshift. 



in all major mergers and in some minor mergers (see Section 2.3). 
Which type of merger dominates the triggering of SMGs? We an- 
swer this question in Fig.|5] where we plot the distribution of stellar 
masses of z = descendants of bright SMGs at z > 1 separated 
into those triggered by minor mergers (blue) and produced by ma- 
jor mergers (red). In this plot, we have weighted the contributions 
to the distribution from the different redshift intervals by the co- 
moving volume per solid angle, so that the final distribution (which 
is normalized to unit area under the histogram) corresponds to what 
would be seen in a flux-limited sample of SMGs over a fixed solid 
angle. We find that 77% are produced by minor mergers, 22% are 
produced by major mergers, and 0.7% are quiescent galaxies (i.e. 
not ongoing starbursts). Minor and major mergers are responsible 
for similar numbers of SMGs for the highest decendant masses, 
but for all other descendant masses, minor mergers predominate. 
These proportions are qualitatively consistent with the result shown 
in Fig. |4] that the bulge-to- total stellar mass ratios B/T of SMGs 
are typically intermediate between pure disk and pure bulge. 



5 PROPERTIES OF SMG DESCENDANTS 

Most of the SMGs detected in flux-limited surveys lie at redshifts 
z=l-4. SMGs at redshift z < 1 account for only a small fraction of 
the population in a flux limited sample due to the relatively small 
volume at low redshift. Model SMGs at low redshifts z < 1 are 
mostly objects with modest star formation rates which have bright 
sub-mm fluxes simply due to their proximity to us compared with 
the bulk of the SMG population, rather than because an intense 
episode of star formation has taken place. In this section, we will 
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Figure 5. The distribution of stellar masses of the descendants of bright 
(Sv > 5mjy) SMGs at z > 1, showing the separate contributions of 
SMGs produced by major mergers (red), minor mergers (blue) and quies- 
cent galaxies (green). The contributions to this distribution from different 
redshifts have been weighted so as to correspond to what would be seen 
in a flux-limited survey over a fixed solid angle. The total distribution is 
normalized to have unit area. 



focus on the descendants of the high-z SMGs, which we define to 
be those with z > 1 on the basis that these have the largest overlap 
with observational samples, and we will investigate the distribution 
in stellar and parent halo mass, morphology and central/satellite 
galaxy classification. 
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Figure 4. Left panels: The predicted distribution of bulge-to-total stellar mass ratio B/T for SMGs selected at different redshifts. Right panels: the B/T 
distribution of the z=0 descendants of SMGs identified at different redshifts. Top panels: bright SMGs selected with Si, > 5mjy. Bottom panels: faint SMGs 
selected with Si, > Imjy. The filled-squares show the median of the distribution at each redshift and the en'orbars indicate the 10% and 90% percentiles. 
The medians are all calculated for volume-limited samples of SMGs at each redshift. 



5.1 Stellar and host halo masses of SMG descendants 

What are the masses of the present-day descendants of SMGs seen 
at different redshifts? In the right panels of Fig. 3, we plot the me- 
dian stellar and host halo masses of the z=0 descendants of the 
SMGs identified at different redshifts. In these panels, the medi- 
ans are calculated for a volume-limited sample of SMG progeni- 
tors, so that descendants are weighted according to the number of 
such progenitors. Note that it is possible for a present-day galaxy to 
have multiple SMG progenitors in different branches of its merger 
tree (see Fig.lTJ. We find that bright SMGs at higher redshift typi- 
cally evolve into somewhat more massive galaxies hosted in more 
massive haloes compared to bright SMGs found at lower redshifts. 
The model predicts that galaxies seen as bright SMGs at the me- 
dian of the observed redshift distribution (z ~ 2) increase their 
stellar mass by one order of magnitude by redshift z=0 ( from 
- 2 X lO^°/i"^A/0 to ~ 2 X lO"/i"^M0). 

In Fig. [6] we show the distributions of stellar and halo mass 
for the z=0 descendants of galaxies which had SMG progenitors 
at 2 > 1. In this figure, and also in Figs. |7] and [S] the distribu- 
tions are calculated for volume-limited samples of SMG descen- 
dants, with each descendant given the same weight whether it had 
one or more SMG progenitors. The top left panel of Fig. [6] shows 
the distribution of stellar masses of the descendants of all bright 
SMGs with z > 1. We separate the galaxies into bulge-dominated 
(B/T > 0.5) and disk-dominated (B/T < 0.5) samples. The 
number of bulge-dominated descendants is nearly double the num- 
ber of disk-dominated descendants. Also, the bulge-dominated de- 
scendant galaxies have a higher median mass. For all of the descen- 
dants, the median stellar mass is 1.5 x 10^^ Mq. The distribu- 



tion of host dark matter halo masses of the z=0 descendants, plot- 
ted in the top right panel of Fig.|6l shows a higher median mass for 
bulge-dominated than for disk-dominated descendants. The median 
dark matter halo mass for all z = descendants of bright z > 1 
SMGs is 2.9 X Mq. For comparison, we also show the 

results for the faint SMGs selected with Si, > ImJy. The median 
stellar and dark matter halo masses of the descendants are about 4 
times smaller than for the bright SMGs, and the distributions are 
wider. In this case, the proportions of disk-dominated and bulge- 
dominated descendants are reversed compared to the bright SMGs, 
with more than twice as many disk-dominated as bulge-dominated 
descendants. 



In Fig.l?] we repeat the same analysis but separate the descen- 
dants into central and satellite galaxies in the halo where they live 
at the present-day. We can see that nearly 80% of the z = de- 
scendants of bright SMGs are central galaxies, while this fraction 
is reduced to 60% for the faint SMGs. We note that satellite galax- 
ies that are descendants of bright SMGs are found only in the more 
massive haloes today. 



We note that lHainline et al.l j2010l) arrived at a similar estimate 
to our own for the masses of the present-day descendants of bright 
SMGs, but starting from observations of SMGs at high redshift, 
and making different assumptions from ours about their subsequent 
evolution down to z = 0. 
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Figure 6. Stellar mass distribution (left panels) and halo mass distribution (right panels) for 2 = descendants of bright z > 1 SMGs (Su > 5mjy, top 
panels) and for descendants of faint z > 1 SMGs (Sy > Imjy, bottom panels). The distributions are for volume-limited samples of descendants, and weight 
each descendant equally whether it had one or more SMG progenitors. The red distribution is for bulge-dominated (B/T > 0.5) descendant galaxies and the 
blue for disk-dominated (B/T < 0.5) descendant galaxies. The black lines represent the sum of the red and blue distributions. The mass distributions are 
normalized to unit area for the full sample. The vertical black arrows indicate the median of the total distribution in each panel. The percentages in each panel 
refer to the percentage of bulge-dominated (in red) and disk-dominated (in blue) galaxies. 



5.2 Morphology of descendants of SMGs compared to other 
similar mass galaxies 

The top right panel of Fig. |4] shows us that the model predicts that 
bright SMGs identified at different redshifts typically evolve into 
highly bulge-dominated galaxies at redshift z=0. According to the 
bottom right panel of Fig. [4] the 2 = descendants of faint SMGs 
at low and intermediate redshifts have intermediate values of B/T. 
Only the higher redshift samples evolve predominantly into highly 
bulge-dominated systems at the present-day. 

Fig. |6] shows that the majority of the descendants of bright 
z > 1 SMGs are bulge-dominated systems, while for faint SMGs 
most descendants are disk-dominated. In Fig.[8]we show (in black) 
the probability distribution of B/T for the z — descendant galax- 
ies of z > 1 SMGs (with each descendant given equal weight 
whether it had one or more SMG progenitors). The top panel shows 
the descendants of bright SMGs, of which 63% of the descendants 
are bulge-dominated (B/T > 0.5) and 37% are disk-dominated 
(B/T < 0.5). Among these, 42% of the descendants are pure 
bulge systems (B/T > 0.9) and only 6% are pure disk systems 
(B/T < 0.1). The bottom panel of Fig. [8] shows the probabil- 
ity distribution of B/T for the descendants of faint SMGs. In this 
case, only 30% of the descendants are bulge-dominated systems, 
with 11% being pure bulge systems, and 12% being pure disk sys- 
tems. 

However, are these results a special characteristic of SMG 
descendants, or are these morphologies simply typical of 2; = 
galaxies with similar stellar mass? To answer this, we select com- 



parison samples (for bright and faint SMGs) composed of galaxies 
at z = with similar stellar mass distributions to the bright and 
faint SMG descendants respectively. In Fig. [8] we plot the proba- 
bility distribution of B/T for the SMG descendants in black, and 
for the comparison samples in red. The descendant and comparison 
samples for bright and faint SMGs show very similar distributions. 
Despite the fact that most of the descendants of bright SMGs are 
bulge-dominated galaxies, their morphologies are not significantly 
different from other z — galaxies with similar stellar masses. 

5.3 Metallicity of descendants of SMGs compared to other 
galaxies 

In the top panel of Fig. |9] we compare the median stellar metal- 
licity (weighted by stellar mass) as a function of stellar mass for 
galaxies which are descendants of bright SMGs (red curves) to that 
of other galaxies of the same stellar mass which are not descen- 
dants of bright SMGs (blue curves). The medians are calculated 
for volume-limited samples of present-day galaxies. The dark and 
light grey shaded regions respectively show the 25-75% and 10- 
90% ranges for the stellar masses of descendants of bright SMGs. 
We see that high stellar mass descendants of bright SMGs have very 
similar metallicities to other present-day galaxies of the same mass, 
while low stellar mass descendants typically have larger metallici- 
ties than high-mass descendants or other galaxies of similar masses. 
For descendants of the median stellar mass, the median metallicity 
is only 20% higher than for non-descendants of the same mass, but 
for the 10% least massive descendants, the difference in median 



© 2009 RAS, MNRAS 000.[Tlfl4l 



10 J. E. Gonzalez et al. 



Stellor moss 



Halo moss 




^ 1-0 

o 

? 0.8 

SI 

^ 0.6 



C7> 
O 



Q. 



0.4 



0.2 



0.0 



c 1 centrol 

S..>5.0mjy gQteHite 



76% 
24% 




10 11 12 
log(M./h-'Me,) 



13 



11 12 13 14 15 
iog(M,„,yh-'M0) 



10 11 12 
log(M./h-'MQ) 



12 13 14 15 

iog(M,„,yh-'M0) 



16 




16 



Figure 7. Similar to Fig|6]but separating the descendant galaxies into central (red) and satellite (blue) galaxies. In this case, the percentages in each panel refer 
to the percentage of central (upper value) and satellite (lower value) galaxies. 



metallicity is a factor 2 or more. It thus appears that for descen- 
dants below the median stellar mass, the starburst which powered 
the SMG was also responsible for chemically enriching the descen- 
dant galaxy above the typical level for galaxies of that present-day 
mass. 

The lower two panels of Fig. |9] make similar comparisons of 
the bulge-to-total stellar mass ratio, B /T, and the halo mass for 
descendants and non-descendants of bright SMGs. We see that the 
most massive 90% of descendants have very similar median B /T 
to non-descendants of the same mass, while the least massive 10% 
have significantly larger B/T than the non-descendants (up to a 
factor 2-3). The halo masses for the most massive 75% of descen- 
dants are very similar to non-descendants of the same mass, but for 
the least massive descendants, the median halo masses are much 
larger (up to a factor ~ 10^). 



6 SMG PROGENITORS OF PRESENT-DAY GALAXIES 

In the previous section we studied the properties of the descen- 
dants of SMGs. We identified SMGs in the model at a particular 
redshift and we found their z=0 descendants. In this section we 
want to study the connection of SMGs to present day galaxies from 
the opposite perspective: given a galaxy today with a particular 
mass, what is the probability that this galaxy had SMG progeni- 
tors? And for dark matter haloes at present-day, what is the proba- 
bility that a progenitor halo hosted at least one SMG? Fig. [TO] an- 
swers these questions, for the 2 = descendants of both bright 
and faint SMGs at z > 1 (top and bottom panels respectively). The 
left panels show this probability as a function of stellar mass, and 
the right panels as a function of halo mass. We see that the prob- 
ability that a galaxy or halo had one or more SMG progenitors at 



z > 1 brighter than a given 850/im flux increases with increas- 
ing mass. From the top left panel, we see that present-day galax- 
ies have a probability of 10% to be descendants of bright SMGs 
{Sv > 5mJy) at a stellar mass of 10"/i-^MQ, and 50% at a 
mass of 6 X 10^^h~^ Mq. For comparison, we see from the lower- 
left panel that 50% of galaxies with present day stellar masses of 
2 X lQ^°h-'^MQ are descendants of faint SMGs (S^ > ImJy). 
We find in both cases (bright and faint SMGs) that at intermediate 
masses (~ 10" -10"/i"^MQ) there are more disk-dominated de- 
scendants of SMGs than bulge-dominated descendants. For haloes, 
we see from the right panels that more than 50% of the haloes 
with masses larger than \Q^'^h~^ Mq are predicted to have bright 
SMG progenitors, and 50% of the haloes with masses larger than 
2 X 10^^/i~^Mq are predicted to have faint SMG progenitors. 



7 CONTRIBUTION OF SMGS TO THE PRESENT-DAY 
STELLAR MASS 

We saw in the previous section that most massive galaxies today are 
predicted to be descendants of SMGs. We found also in Section 4. 1 
that the median of the SMG phase duration is 0. 1 IGyr. Since in the 
literature a high star formation rate (SFR) is traditionally inferred 
for the bright SMGs (based on the assumption of a normal IMF), 
we want to know what contribution to the present-day stellar mass 
is actually produced in the SMG phase. 

In Fig. [TT]we plot the predicted evolution of the cosinic SFR 
density. We show both the total, as well as the contributions from 
bursts and quiescent disks, and the contributions from galaxies in 
a bright {Si, > 5mJy) or faint (5*1, > ImJy) SMG phase. We see 
that the contribution to star formation from quiescent disks domi- 
nates at redshift z < 3.5, but bursts dominate at higher redshifts. 
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Figure 10. Left panels: fi'action of present-day galaxies by stellar mass that are predicted to be descendants of bright SMGs (Sv > 5mjy, top panel) and 
descendants of faint SMGs (Sv > Imjy, bottom panel). Right panels: fraction of present-day halos by mass that are predicted to have halo progenitors 
that hosted at least one blight SMG (top panel) or at least one faint SMG (bottom panel). These are computed for z > 1 SMG progenitors. Left panels: 
bulge-dominated (B/T > 0.5) descendants in red and disk-dominated (B/T < 0.5) descendants in blue. Right panels: bulge-dominated (B/T > 0.5) halo 
central galaxy in red and disk-dominated (B/T > 0.5) halo central galaxy in blue. The black line represents all of the SMG descendants (the sum of the blue 
and red lines). 



as found earlier by 'Baug h et al] ( l2005h . At redshifts z = 2 — A, 
around 1% of the cosmic star formation is produced by bursts in 
bright SMGs, and around 10% by bursts in faint SMGs. The quies- 
cent SFR density is typically one order of magnitude smaller than 
the burst SFR density for bright SMGs at high redshift. For faint 
SMGs, the burst SFR density dominates at redshift 2 > 2, but qui- 
escent star formation becomes more important at lower redshift. 
At low enough redshift {z <g 1), star formation at any level in 
a galaxy is sufficient to produce a sub-mm flux above the bright 
or faint thresholds, which is why the quiescent and burst contribu- 
tions from bright and faint SMGs converge on the total quiescent 
and burst contributions &l z = Q. 

In the top panel of Fig.[T2]we plot the comoving stellar mass 
density as a function of redshift, showing both the total, and the 
contributions to this from star formation in the quiescent and burst 
modes, and the contributions from star formation during bright and 
faint SMG phases at z > 1. Note that the stellar mass densities 
we plot are the values after allowing for recycling of gas to the 
ISM from dying stars, using the recycled fractions for the two IMFs 
given in Sec. |2.3l (recall that the recycled fraction for the top heavy 
IMF is close to unity). The lower panel shows the different contri- 
butions to the stellar mass density plotted as fractions of the total 
stellar mass density at that redshift. We see from this that overall, 
bursts of star formation contribute a total of 5% to the present-day 
stellar mass density. Star formation (burst-l-quiescent) in the faint 
SMG phase at 2 > 1 contributes 2% to the present-day stellar mass 
density, while the bright SMG phase aX z > 1 contributes only 
0.06%. The contribution of bursts and SMGs to metal production 



are however much larger, due to both the larger yield and larger re- 
cycled fraction for the burst IMF compared to the quiescent IMF. 
The model predicts that all bursts over the history of the universe 
produced 70% of the metals existing today, while bright SMGs at 
2 > 1 produced 0.8% of the metals existing today. 

Finally, we examine in more detail the contribution of the 
mass produced by all bursts and by bursts in the bright SMG phase 
as a function of stellar mass. In Fig. [T3] we plot the fraction of 
the present stellar mass (allowing for recycling) produced by all 
bursts in all galaxies (red) and only in descendants of bright SMGs 
(black). We also plot the contribution to the stellar mass produced 
by bursts and quiescent star formation in the bright SMG phase 
for all present-day galaxies (green) and only for the descendants of 
bright SMGs (blue). We see from this plot that on average bursts 
contribute 3-8% to the present-day stellar masses of galaxies over 
the whole range ^ Iff ~ lO^^hT^ Mq. If we look only at galax- 
ies which are descendants of bright SMGs, then we find an average 
fraction that is the same at the highest stellar masses, but increases 
with decreasing mass, reaching 60% for the least massive descen- 
dants (~ 5x 1O^^~^M0). However, the total stellar mass produced 
in the bright SMG phase at 2 > 1 is less than this. For descendants 
of bright SMGs, the fraction varies from 30-40% for the least mas- 
sive descendants (~ 5 x lO'-^ft^^M©), to 0.2-0.4% for the most 
massive ones (~ IQ^^h-^MQ). The average contribution to the 
current stellar mass of all galaxies from bursts in the bright SMG 
phase increases with stellar mass, but is always below 0.5%. For a 
present-day stellar mass of L5 x lO^^/i'^A/©, which is the median 
descendant mass of bright SMGs (see Sec. 15. It , the contribution of 
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Figure 8. Probability distribution of the bulge-to-total stellar mass ratio 
(B/T) for z = Q descendants of bright z > 1 SMGs (top panel, black) 
and for descendants of faint z > 1 SMGs (bottom panel, black). The dis- 
tributions are for volume-limited samples of descendants, and give each 
descendant the same weight whether it had one or more SMG progenitors. 
The dashed red lines show the B/T distributions of the comparison sam- 
ples formed by galaxies at z = with a similar stellar mass distribution to 
the bright and faint SMG descendants. All distributions are normalized to 
unit area under the histogram. 



Stellar mass produced in the bright SMG phase to the present stellar 
mass is only 1-2% even for descendants of bright SMGs, and falls 
to only 0.3% if we look at all galaxies of that mass. 



8 CONCLUSIONS 

We have investigated predictions for the properties and descendants 
of submillimetre galaxies (SMGs ) in a theoretical m odel of galaxy 
formation in a ACDM universe jBaugh et alj|200^ . This model, 
which incorporates a top-heavy IMF in starbursts, has previously 
been shown to reproduce the observed number counts and redshifts 
of SMGs, as well as a wide variety o f other galaxy properties at low 
and high redshift. We emphasize that lBaugh et al J found that the in- 
clusion of a top-heavy IMF in starbursts was essential in order to 
match both the numbers and redshifts of SMGs and the present-day 
galaxy luminosity functions in a single model. If instead the IMF 
was assumed to be universal and of solar neighbourhood form, then 
by adjusting parameters for feedback and other processes, either the 
present-day galaxy luminosity functions could be reproduced, but 
the sub-mm counts underpredicted, or the sub-mm counts could be 
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Figure 9. The three panels show the median stellar metalhcity, bulge-to- 
total stellar mass ratio, and halo mass vs stellar mass for present-day de- 
scendants (red hne) and non-descendants (blue line) of bright SMGs. The 
medians are calculated for volume-limited samples of present-day galaxies. 
The error bars show the 10-90% range in bins of stellar mass. The dark grey 
shaded region shows the 25-75% range for the stellar masses of descendants 
of bright SMGs, while the light grey shading shows the 10-90% range. 



reproduced, but at the price of predicting far too many high lumi- 
nosity (and high stellar mass) galaxies at the present day. 

We find that bright SMGs in the model, which we define to be 
objects having 850/im fluxes brighter than 5 mjy at some redshift 
z > 1, are typically the progenitors of massive galaxies today, with 
a median stellar mass ~ 2 x lQ"h~^ Mq for the descendants. 
The present-day descendants of these bright SMGs are predicted 
to be mainly bulge-dominated systems, but the bright SMG phase 
does not appear to play a dominant role in building the bulge, since 
their morphologies are similar to those of other present-day galax- 
ies with similar stellar masses. The descendants of bright SMGs are 
also predicted to live in the most massive haloes today. Only the 
most massive objects in the present-day universe are predicted to 
be predominantly descendants of bright SMGs - we find that more 
than 50% of galaxies with masses above 6 x 10"/i"^Mq andha- 
los with masses above 10^^/i~^Mq have one of more progenitors 
which hosted bright SMGs at some stage in their evolution. 

The model also predicts that the typical stellar masses of 
bright SMGs at 2: = 2 are ~ 2 x lO^°fc~^M0 similar to the new 
observational estimates bv lHainlineetal.l ( l201(]|) . and that they are 
typically mildly bulge-dominated in stellar mass, consistent with 
the observations by Swinbank et al. (2010). The typical duration of 
the bright SMG phase is predicted to be ~ O.lGyr, again consis- 
tent with observational estimates. 

In our model, starbursts account for around 30% of all star 
formation over the history of the universe. However, stars formed 
in galaxies which are detectable as bright SMGs at 2; > 1 are pre- 
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Figure 11. Evolution of the cosmic star formation rate per unit comoving 
volume. The black line shows the total total star formation rate, while the 
green and blue lines show the separate contributions from bursts and quies- 
cent star formation respectively. The red and orange lines respectively show 
the star formation rate in bright and faint SMGs, for which the separate 
contributions from bursts and quiescent star formation are shown by dashed 
and dotted lines, respectively. 



dieted to aecount for only 0.06% of the stellar mass density (in 
live stars and stellar remnants) today. There are two reasons why 
the latter fraction is so small. Firstly, the recycled fraction for the 
top-heavy IMF (i.e. the fraction of the initial stellar mass which 
is returned as gas to the ISM when stars die) is > 90%, much 
larger than for a solar neighbourhood IMF, so that less than 10% of 
the stellar mass formed in bursts remains in stars and remnants to- 
day. Secondly, only a small fraction of ongoing starbursts at 2 > 1 
have high enough star formation rates and luminosities to be de- 
tectable as bright SMGs. For present-day galaxies which are de- 
scended from bright SMGs at 2; > 1, the fraction of the current 
stellar mass which was formed in the bright SMG phase is larger, 
ranging from ~ 30% in low mass descendants to ~ 0.3% in the 
highest mass descendants. However, for median mass descendants, 
the mean fraction formed in the bright SMG phase is still only 
-3%. 

In summary, our model predicts that SMGs are the progeni- 
tors of massive galaxies today. However, most of the stellar mass in 
these systems is built up by quiescent star formation and then as- 
sembled in galaxy mergers, making the contribution of long-lived 
stars formed during the SMG phase typically very small. The latter 
is contrary to the view commonly expressed in the observational lit- 
erature on SMGs. Thus in our model, the SMG phase mainly serves 
as a "beacon" for gas-rich mergers in the progenitors of massive 
galaxies at high-z. 



ACKNOWLEDGEMENTS 

We thank the referee for their constructive report and suggestions 
which helped to improve the paper. JEG acknowledges receipt of a 
fellowship funded by the European Commission's Framework Pro- 
gramme 6, through the Marie Curie Early Stage Training project 
MEST-CT-2005-021074. This work was supported in part by the 



10 



10 



u 

CL 



10' 



10' 



10' 




N 



^ 10"' 

N 

1 10-' 

10"^ 

10-^ 
10"^ 



N 



quiescent 
burst 
foint SMGs 
bright SMGs 



|| I I I I 



1 



2 3 4 
z 



Figure 12. Top panel: comoving stellar mass density as a function of red- 
shift. The black line shows the total, while the blue and green lines show the 
contributions to this from quiescent and burst stai' formation respectively, 
and the red and orange lines show the contributions from star formation in 
bright or faint SMG phases. For the SMGs, we only show the contributions 
from SMG phases at 2 > 1. Bottom panel: similar to the top, except that 
the different contributions to the stellar mass density are plotted as fractions 
of the total stellar mass density at that redshift. 
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